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I. INTRODUCTION 

The wide-spread interest in quantum information process- 
ing in recent years has been a critical driving force in the re- 
search of electron spins localized in semiconductor quantum 
dots (QDs).' While these two-level systems are attractive can- 
didates for implementation of scalable systems due to their 
compatibility with conventional microelectronic technology, 
their quantum control at the single- and few-qubit level re- 
mains a significant experimental challenge. Furthermore, as 
all solid state systems, they are inherently less isolated from 
their environment as compared with atomic systems. 

An attractive platform to study quantum control and the re- 
lated problem of decoherence is the system of gate-defined 
lateral QDs, in which several of the major breakthroughs in 
spin qubit technology have emerged in recent years. While 
isolating a single electron in a QD was achieved only in 2000,^ 
rapid progress has been made since then. Long singlet-triplet 
relaxation times of the order of milliseconds were measured 
for a single dot,^ and a lower bound on the spin coherence 
time (dominated by pure dephasing) exceeding 1/is was es- 
tablished, using spin-echo techniques in a double dot systemii 
The single-spin relaxation and decoherence time scales have 
since been pushed to the order of 1 s^ and 0. 1 ms'' respectively. 

The relative isolation of QD electron spins, which is in- 
dicated by these long coherence times, renders their ma- 
nipulation and readout particularly challenging. This is ac- 
complished by using Pauli spin blockade to convert spin to 
charge information so that fast measurement of spin states 
becomes possible. In addition, coherent exchange of two- 
electron spins in a double dot system^, and driven Rabi os- 
cillations of single electron spins using oscillating magnetic"^ 
and electric" fields have been demonstrated as well. 

Electron spin relaxation via spin-orbit interaction was 
shown to be an insignificant decoherence channel,'^ and it 
has been generally accepted that the nuclear spins in the 
surrounding host material are the main source for the elec- 
tron spin decoherence in III-V host materials such as GaAs 



and InAs»i^ This has led to intensive experimenta l'^i i^^ii^i^^-'"^-'^ 
and theoretical ' ^,20,21 .22,23,24.25,26 studies of the nuclear en- 
vironment, and various proposals for alleviating its ad- 
verse effects on the electron spin qubit, among which dy- 
namical nuclear spin polarization was suggested2i2i22 and 
demonstrated J ' 

In contrast, the effects of the charge environment on QD 
spin qubits have only recently started to receive some theo- 
retical attention. '^'^^'^"^ Charge noise in lateral gated devices 
can originate from various sources. Suggested mechanisms 
include gate leakage currents via localized states, charge traps 
near the quantum point contacts (QPCs), donor centers near 
the gate surface, Johnson noise from the gate electrodes, and 
switching events in the doping layer, typically located at an 
interface 100 nm below the surface. "^'^^ Measurement of the 
background charge fluctuation in GaAs quantum dots has 
shown a linear temperature dependence characteristic of 1// 
noise. Random telegraph noise in GaAs lateral gated struc- 
tures was measured and characterized by Pioro-Ladriere et 
al. ^ This noise was attributed to electrons that tunnel from 
the gate and are trapped near the QPC, causing fluctuations in 
the conductance with typical frequency of 1 Hz. Applying 
a positive gate bias during the device cooldown significantly 
reduces the noise by reducing the density of ionized donors 
near the surface, thereby suppressing the electron tunneling. ^2, 
Furthermore, background charge fluctuations were suggested 
as a possible source for the bistable behavior observed in the 
coupled electron-nuclear spin system,'^ and telegraph noise 
induced by the QPCs was also measured recently in double 
and triple coupled QDs.^^ 

Generally, single-spin qubits in solids rely on the exchange 
interaction to perform fast two-qubit operations. Furthermore, 
a number of recent works have utilized two-spin singlet and 
unpolarized triplet states in biased configuration to encode a 
logical qubit, which offer better control as compared with sin- 
gle spin states.^** However, such exchange-coupled spin qubits 
are vulnerable to dephasing induced by charge noise, since 
exchange coupling is electrostatic in nature, and singlet and 
triplet states generally have different charge distributions.— 



2 



In the same spirit, the effects of charge noise on the coher- 
ence of spin qubits in Silicon double dots were studied very 
recently.— In addition, the effects of a single, randomly po- 
sitioned, charge impurity on a three-spin encoded qubit in a 
triple QD were studied by calculating the impurity-induced 
changes in the qubit orbital levels.— Finally, electron-phonon 
interaction can also lead to dephasing in an exchange coupled 
double quantum dot because two-spin singlet and triplet states 
have different charge distributions 

There are many types of charge impurities and defects that 
can generate electrical fluctuations that affect spin qubits in 
solid states. In this paper we carry out a microscopic cal- 
culation focused on the Coulomb coupling between a biased 
two-spin qubit and a nearby trapped charge fluctuator repre- 
sented by a two-centered two-level-system (TLS), utilizing 
a multipole expansion up to and including the quadrupole- 
quadrupole order One scenario for such a two-center de- 
fect may be for an electron to be trapped around two donor 
nuclei that have potential wells somewhat lower than other 
donor nuclei nearby, so that this electron would oscillate be- 
tween the sites until the charge motion is relaxed by the back- 
ground charge fluctuations or phonon emissions. Using a mas- 
ter equation formalism we use the calculated qubit-TLS cou- 
plings to study the dynamics of the open system that is formed 
by the spontaneous emission of the TLS coupled to a reservoir 
Thus we obtain quantitative estimates of the decoherence and 
dephasing effects on the spin qubit during various gate opera- 
tions, and when idle. This analysis enables us to determine op- 
timal working points at which the qubit' s sensitivity to charge 
fluctuations is reduced. 

It is important to note that this work is only an initial step in 
the quantitative analysis of the effects of charge fluctuations 
on spin qubits. The focus here is on a quantitative evalua- 
tion of the qubit-TLS entangling term. We are particularly in- 
terested in clarifying how the TLS-qubit coupling could lead 
to qubit decoherence due to the background charge fluctua- 
tions, with the TLS acting as an intermediary between the spin 
qubit and the charge environment. The TLS coupling to the 
environment is dealt at a rudimentary level, serving only to 
demonstrate the applicability of the presented theory in es- 
timating charge-induced spin decoherence. Building on the 
results given in this paper, the next step should benefit from 
the extensive work that has been carried out in recent years on 
charge-environment-induced decoherence in superconducting 

qubitS.^^-^3,44,45,46,47,48 

The paper is organized as follows. In section |II] we de- 
rive the Coulomb coupling between the qubit and the TLS, us- 
ing the Hund-Mulliken approach to calculate the qubit orbital 
states and a multipole expansion for the Coulomb interaction. 
In section |III] we use these results to study the decoherence 
effects due to charge fluctuations mediated by the qubit-TLS 
coupling. After deriving the master equations for the system 
density matrix in section llll Al we present and analyze in sec- 
tions |IirB]|lirc] the resulting dynamics during various single- 
qubit operations for singlet- triplet qubits. In section HIlDl we 
discuss a convenient working point at which the effective ex- 
change energy is zero and quantify the dephasing time. A 
summary of our results and a brief discussion on possible ex- 



tensions of this work are given in section |IV] In Appendix A 
we calculate the system concurrence, showing the conditions 
for qubit-TLS entanglement. In appendix B we provide de- 
tails of the qubit's orbital Hamiltonian. Appendix C lists the 
full expressions for the qubit-TLS coupling terms, and finally 
appendix D presents an analytical solution to the master equa- 
tion for the case of no TLS tunneling. 

n. QUBIT-TLS COULOMB COUPLING 

We consider the Coulomb interaction between a qubit 
formed from the singlet and unpolarized triplet spin states 
of two electrons in a double dot and a nearby two level sys- 
tem (TLS), assuming no qubit-TLS tunnel coupling. To prop- 
erly describe a biased double quantum dot, we use the Hund- 
Mulliken model to calculate the qubit orbital states. With 
Coulomb interaction being spin-independent, the interaction 
Hamiltonian can generally be written in the form:"^** 

Hint ^ -aaj (g) I'^ - (31'^ (g) af + jaj (g) af , (1) 

where we have 

a = - (y™ + - _ ySL^ 
4 ^ ^ 

/3 = i (y™ _ ySR ^ yTL _ ^SL^ 

Here are the Coulomb matrix elements, where the left su- 
perscript denotes the qubit state (Singlet/Triplet) and the right 
one denotes the TLS state (Left/Right). While the a coupling 
should not directly affect the qubit spin state, the /3 coupling 
effectively renormalizes the qubit exchange energy. The 7 
coupling acts to entangle the qubit and the TLS and therefore 
leads to qubit spin decoherence when the TLS is coupled to 
a larger reservoir representing the background charge fluctua- 
tions. Appendix A formalizes this last statement, showing that 
the concurrence^" of qubit-TLS system under the time evolu- 
tion of Hint is nonzero only for a nonzero 7 coupling. 

The TLS's we consider are sufficiently removed from the 
spin qubit so that there is no exchange coupling between the 
spins and the single electron in the TLS. With no electrons 
tunneling between the qubit and the TLS, the two charge dis- 
tributions are separated in space, and the Coulomb interaction 
between them can be described systematically using a multi- 
pole expansion approach In the following we use this model 
to evaluate the Coulomb coupling terms a, f3, 7. 

A. Hund-Mulliken approach for the qubit orbital Hamiltonian 

The sensitivity of an exchange-coupled spin qubit to a re- 
mote charge fluctuator comes from the different charge distri- 
butions the singlet and triplet states have.''' Thus we first con- 
struct the two-electron orbital states by extending the Hund- 
Mulliken approach^^ to a biased dot configuration. 
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We start by approximating the orbitals for the two quantum 
dots by those of two harmonic wells centered at ±ax 



jmuj ±i- 



X{z), (3) 



where uj = huQ. Here loq is the harmonic oscillator frequency, 
and the magnetic compression factor, 6, is given by 6 = 
a/1 + cj|^/ijJq with the Larmor frequency lol — eB/2mc. 
The phase factor in Eq. (O involving the magnetic length 
Ib = \/ hc/eB results from a gauge transformation, and 
Aa; = eE/mujQ is the orbital shift due to the electric field. 
The z direction wavefunction is taken as the ground state (with 
associated energy Eqz) of a finite potential well Vz of width 



COs(fceZ) 

fc^ ^-k(\z\-L,/2) 



< Lz/2 
> Lz/2 



(4) 



with fee = ^{2m,/h^)Eoz, n = ^i2m,/h^)iVz- Eoz), 
Nz = {Lz/2 + l/n)-^/\ 

To simplify the Hund-Mulliken calculation, the single par- 
ticle single-dot states are orthonormalized, ■4>±a = N'{(f>±a — 
g(l)^a), where g = (1 - Vl - s^)/s, s = {(f>a\(f>-a) 
g-d (2fc-i/&) jg jjjg wavefunction overlap (d ~ a/as, with 
a B the Bohr radius associated with the harmonic QD confine- 
ment potential), and Af = {1 — 2sg + The orthonor- 
malized orbitals are then used to construct 4 two-particle 
states: the two doubly occupied singlets, S'(2, 0) = V'-aV'-a, 
S{0,2) — ipa'tpa, the separated singlet state, 5'(1, 1) = 
{tp-ai^a + ipa'ip-a)/ V^, and the separated triplet state, T = 
r(l,l) = [ip^ai'a - i-'ai'^a)/\/^- (wc ncglcct the doubly 
occupied triplet states as their energy is typically much higher 
for the gate-defined structure we have in mind.^''"'''*) In the 
basis of these two-particle states the orbital Hamiltonian is 
given as 
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(5) 



where the diagonal elements include the Coulomb interac- 
tions, and tn and X are the single- and double-hopping ma- 
trix elements, respectively. Calculational details of the or- 
bital Hamiltonian in a biased configuration are given in ap- 
pendix B. Figure[T]shows the energy diagram near the (1 , 1) to 
(0, 2) charge transition, which results from diagonalization of 
Eq. (|5]l, where we also included the polarized triplet states T± 
splitted by the Zeeman interaction, Hz = guB^^- Si=L r 
with g = —0.44 and /.is the Bohr magneton. For this figure 
and throughout this work we have considered B — 100 mT 
(Ez = 2.5/ieV), dot confinement luq — 3meV (as — 19.5 
nm), and half interdot distance (in gb units) d = 2.8, cor- 
responding to the experimental parameters in typical gate- 
defined double dot systems."*'"' The bias shift in the figure is 
normahzed to the Bohr radius: 
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FIG. 1: (Color online) Orbital energy diagram for the double dot 
near the (l,l)-(0,2) transition vs the dimensionless bias shift Ax. 
Shown are the Hybridized singlet states (black curves) and split (1,1) 
triplet states T- (dash-dotted green), To (solid blue), and T+ (dashed 
red). The inset shows the exchange energy, calculated as the energy 
difference between the separated triplet state and the lowest lying 
hybridized singlet state. Here and throughout the paper B = 100 
mT, d = 2.8, and a;o = 3meV 



and it is proportional to the interdot bias gate potential. 

Since the orbital Hamiltonian, Eq. (|5]l, does not connect the 
triplet state with any of the singlet states, the combined two- 
particle orbital-spin triplet state can be written as: 



Ti + iT) 

^/2 ■ 



(6) 



The diagonalization of the singlet 3x3 block of iJoib yields 
a hybridized singlet state that is predominantly the separated 
singlet, ^(Ijl) at negative or zero bias.^' In the basis of 
the three singlet states {5'(2, 0), S{0, 2), S'(l, 1)} the lowest- 
energy orbital-spin singlet state is 



I Ti - it) 




= TVs [ailjj^alp-a + a2'>Pa^a 

I Ti - it) 



V2 



(7) 



where Ms = 1/ \Jl + a\+ a\, and ai , 02 are the 
5(2,0),S'(0,2) components of the lowest lying singlet eigen- 
state of the orbital Hamiltonian Eq. (|5]l. The exchange energy 
is defined as the difference between the triplet and this singlet 
state (see inset of Figure[T]|. 

The two-electron states can be expressed in terms of the 
single particle orthonormal states -01 = (V'-a + ^a)/>/2, 
■02 = (V'-a — ^a)/ \/2, which are more convenient when cal- 
culating the qubit-TLS couplings. The two-particle orbital- 
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spin states are given by the slater determinants: 
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V2 

I Ti - it) 
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= [^^1(1)^2(2)1 U) - V'2(l)^i(2)UT)] , 

where the number index on the left-hand-side denotes the first 
or second electron. The triplet and hybridized singlet states 
can be built as 

\To) = ^ (V'-aV'a - M-a) ^^^^^ 



-^[|llTl2i> + |lul2T> 



(8) 
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(9) 



Notice that the last term in Eq. (|9]l vanishes for unbiased dou- 
ble dot (ai = 02). In what follows, these states will be used 
to calculate the qubit-TLS coupling terms, Eqs. (|2]i. 



B. Multipole expansion for the qubit-TLS interaction 

We model the TLS as a single electron moving in a double 
well, each of which has a wave function similar to those of the 
qubit orbitals, Eq. Q 



1 



-[{xTaTf+y^ 



^/ttDt 



Xt{z), 



(10) 



are the electron charge density operators for the qubit (Q) and 
TLS (T), where i,j £ {1,2} denote the qubit orbital state 
(symmetric or antisymmetric combination), and fc, I e {L, R} 
denote the TLS state. We use the dielectric constant for GaAs, 
e = 13.1, and consider only static dielectric constant for 
screening, since we assume the space near the double dot is 
completely depleted (i.e. no nearby 2DEG). In addition, we 
take the TLS inter-site distance to be sufficiently large so as 
to have a relatively small tunnel coupling, limiting our study 
to slow TLS's. We can therefore neglect contributions to the 
qubit-TLS coupling coming from off-diagonal TLS charge 
densitiesi^ To reduce clutter we thus write fijk = fijkk- 

The Coulomb matrix elements of interest are V'^'^ = 
{Tk\C\Tk), and F^'^ = {Sk\C\Sk), where T,S are the 
triplet and singlet states given in Eqs. ([8j-(|9]l. Assuming the 
creation operators for the electrons in the QDs commute with 
those in the TLS (i.e., no tunneling between the qubit and the 
TLS) we find 

(Tk\C\Tk) = /life + /22fe 

(^fc|C|5fc) =/iife + /22fe+AA|(ai+a2) (13) 
\/2 (/life - /22fe) + (ai - 02) (/l2fe + /21fc) 



We calculate the Coulomb interaction terms fijk by eval- 
uating the electrostatic energy associated with placing the 
TLS charge distribution in the potential "I>g, that is due to 
the qubit charge distribution, expanding the latter in spherical 
harmonics^^ 



/.,fe= / dr$3(r)p^'=(r) 
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1=0 



9q P4 • r 1 



'Qlm ^5 



(14) 



where Dt is the Bohr radius of the (identical) TLS cen- 
ters, ot is half the distance between them, and xt{z) is the 
TLS ground state z wavefunction similar to that of the qubit, 
Eq. with potential , and width L^.^^ 

The most general Coulomb interaction operator between 
the qubit and the TLS is given by 



fi]ki = / drdv 



e\v - 



(11) 



where e is the dielectric constant, and , Pg , Qg are the 
charge, dipole, and quadrupole electric moments, respec- 
tively, associated with the qubit charge distribution. Combin- 
ing this with the Taylor expansion for the potential 



$(r) = $(0)+r-V$(0) + i^ 



rirr, 



92$ 

dridr„ 



■(0)- 



(15) 



where 



p^4{v) = e4*(r)</,^(r), 



(12) 



and using q^^V^.Q^ to denote the charge, dipole, and 
quadrupole electric moments, respectively, of the TLS charge 
distribution, we obtain /^fe up to and including quadrupole- 
quadrupole order: 
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FIG. 2: qubit-TLS system geometry. 

In Eq. ( fT6] l R = i? (sin6' cos sin6' sin(/<, cos (?) is the vec- 
tor connecting the qubit and TLS centers, and the two dots lie 
along the X-axis. The centers of the TLS are aligned along 
the axis xt = sin^T cos^Ti + sin^T sin^Ty + cos 9tz, so 
the angular dependence of the qubit-TLS interaction is spec- 
ified by the four angles {6, ip, 9t, (t>T)- The system geometry 
is depicted in Figure|2] 

The electrical monopole (charge) for both the qubit and 
TLS is just e.^^ The qubit dipole moments are, by construc- 
tion, in the x direction and are found to be 



Q 



drrp, 



•q 



eaa(lTg)^Aa; 



(l±s)x, (17) 



where the upper (lower) sign corresponds to i = 1 (i = 2). 
For the mixed qubit orbital Coulomb matrix elements we only 
need to consider the sum fi2k + /21/c (see Eq. (T3[ for the 
off-diagonal matrix element). In the case of the qubit dipole 
moments this sum is 



p21 



2eaBd{l - g^) 
1 - 2sg + g2 ■ 



(18) 



Notice that for unbiased dots (Ax — 0) the diagonal qubit 
dipole moments, Pq^ , Pq^ vanish, since in this case the dou- 
ble dot is symmetric ((/)_a(— x) = 4'a{x)) thus the charge dis- 
tribution has mirror symmetry around the y ~ z plane and the 



dipole moment is identically zero. The introduction of bias 
allows one of the double occupied states to mix more strongly 
into the ground singlet state so that a finite dipole moment 
emerges. 

Using the TLS wavefunctions Eq. ( fTOl i we find the TLS 
dipole moment as 



(19) 



where the plus (minus) sign corresponds to k — R{k — L). 
The quadrupole moments Q\l-,-^ — J dr{3rirm ~ r^Sim)p''-' (j^) 
are found to have only diagonal elements. For the qubit charge 
distribution they are 



l~2sg + 52 
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In Eqs. ( |20] | the upper (lower) sign refers to i = 1 (i 

.2,2 



(21) 
2), and 
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Similarly, the quadrupole moments of the TLS charge distri- 
bution are 

Qryy^e (^-a| + ^ - (23) 

where I'^j. = J dzz^Xrlz) is the square extension of the z- 
direction TLS wavefunction, similar to Eq. (l22l l. Note that 
the TLS quadrupole moment matrix elements are given in the 
rotated frame [xt, yr, zt) and have the same values for k — 
L,R. 

Using the above results we obtain the qubit-TLS coupling 
terms, Eqs. (|2|i to quadrupole-quadrupole order where the 
nonvanishing contributions are 

a = Uqd + add + ctQd (24) 

/3 = /3<ig + /3q, + PdQ + PQQ (25) 

7 = 7<id + iQd- (26) 

The first (second) subscript in each term denotes contribution 
from the particular multipole moment: monopole (q), dipole 
(d), and quadrupole (Q) of the qubit (TLS) charge distribu- 
tion. The explicit expressions for the various coupling terms 
are rather lengthy and are deferred to appendix C. Note that 
the dipole-charge, (3dq, dipole-dipole, add, Idd, and dipole- 
quadrupole, fidQ contributions are nonzero only for a biased 
dot configuration when the qubit dipole moment is nonzero. 

C. Qubit-TLS coupling terms 

In order to present graphically the qubit-TLS Coulomb in- 
teraction terms, Eqs. (|24]|-(|26]|. we consider a generic system 
geometry (i.e. no a-priori knowledge of the relative orienta- 
tions of the two subsystems). Since giving statistics of var- 
ious parameters does not shed clear light on the interaction, 
we give a representative value of the qubit-TLS interaction 
by averaging over all possible values of qubit-TLS orientation 
parameters [9, (j), 6t, (t>T)- To obtain the correct total values 
of a, /3, 7, this angular averaging should be performed after 
the addition of the individual contributions from the multipole 
expansion (i.e., the terms given in Eqs. ( |C-lt -( IC-9t ). Figure 
|3] shows the angle-averaged values of the coupling terms vs 
qubit-TLS distance at the bias point corresponding to the sin- 
glet anticrossing (see Figure [TJ. In addition to the double-dot 
parameters given in Fig.[Tl we use here and throughout the pa- 
per, unless specified otherwise, QD thickness = 5 nm, and 
vertical confinement potential Vz = 500 meV. For the TLS, 
we take = 3 nm, = 100 meV, TLS center Bohr radius 
Dt = 5 nm, and half TLS centers distance = 20 nm. The 
latter are chosen to characterize (5-doped dopants in the insu- 
lator with a typical small radius and a fairly large inter-center 
distance. 

The figure demonstrates the convergence of the multipole 
expansion for each of the three couplings, as R increases. For 
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FIG. 3: (Color online) Angle-averaged Coulomb coupling terms vs 
qubit-TLS distance at double-dot bias corresponding to the (1, 1) to 
(0, 2) charge transition (Ax = 0.405). (a) a terms; (b) /3 terms; (c) 
7 terms. 



the above set of parameters we expect higher order contribu- 
tions in the multipole expansion to be insignificant for qubit- 
TLS distances exceeding 100 nm (for a), 40 nm (for (J), and 
30 nm (for 7). As explained below, it is the 7 coupling that 
is responsible for the spin qubit decoherence effects, thus we 
expect our results to be accurate down to i? = 30 nm. 

Figure |4] shows the Coulomb couplings and the qubit ex- 
change energy as functions of the QD bias shift for qubit-TLS 
distances of R = 30, 80, 200 nm. While the bias dependence 
of the a coupling is minimal, both f3 and 7 strongly depend on 
the qubit bias, suggesting that the qubit is substantially more 
susceptible to decoherence due to charge fluctuations at and 
above the anticrossing point where the S{Q, 2) component in- 
creases significantly in the ground singlet state and the qubit 
charge distribution acquires a strong dipole component. It is 
seen that for our parameter choice, 7 becomes comparable to 
the exchange energy at i? < 100 nm, below which we an- 
ticipate sizable spin decoherence effects due to charge cou- 
pling. Inspection of the leading terms in the f3 and 7 couplings 
(Eqs. ( IC-4b -( IC^ ) shows that only the latter scales with the 
TLS centers distance qt, thus the ratio f3/j decreases with 
the charge fluctuator size, leading to increased qubit-TLS dis- 
tances at which the qubit is susceptible to decoherence. We 
note that the qubit-TLS distance at which the Coulomb terms 
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FIG. 4: (Color online) Angle-averaged qubit-TLS coupling terms as 
functions of the double-dot bias for qubit-TLS distances of: (a) R = 
30 nm, (b) R — 80 nm, and (c) R = 200 nm. The qubit exchange 
energy is shown for comparison. All other system parameters are the 
same as those used in Figure [5] 



become appreciable roughly scales linearly with the size of 
the dots, which is consistent with the basic characteristics of a 
multipole expansion. 

The angular averaging procedure that was used to produce 
figures [3] and |4] was tested by randomly taking values for 
{9, ip, 9t, 4>t) and using these random 4-vectors to calculate 
the interaction terms. This calculational mode is useful in later 
evaluation of decoherence effects such as gate errors. We then 
calculate the error (or any other decoherence effect) for many 
randomly selected qubit-TLS geometries and average at the 
end of the calculation. The results shown in Figs. [3] and |4] 
were reproduced to within an error < 1% by averaging over 
10,000 random runs. 



III. EFFECTS OF CHARGE FLUCTUATIONS ON A 
DOUBLE-DOT SPIN QUBIT 



of the charge coupling on dephasing of the spin qubit and on 
the fidelity of specific single-qubit operations. 



A. Master equation for the qubit-TLS system 

We consider the master equation describing the qubit-TLS 
system, with the TLS coupled to a reservoir that results in its 
spontaneous emission 



^ = -l[n,p]+j:[^L,pL]-{L]L,,p} 



(27) 



Here p — pT ® pq is the 4x4 density matrix of the qubit- 
TLS system, Lj are the Lindblad operators, and the qubit-TLS 
Hamiltonian is 



n 



Br 



(28) 



where 7ii„t is given in Eq. ([T]i, and Q (T) superscript de 
notes an operator on the qubit (TLS) subsystem. In Eq 



B 



Q 



■^{Sh, 0, J), with 6h the magnetic field inhomogeneity 
between the dots arising from either application of an inho- 
mogeneous B, different g factors in the two dots, or inhomo- 
geneous nuclear polarizations. By = {tr, 0, ojt) where ujt is 
the TLS level splitting and tr is the tunnel coupling between 
the two centers. The latter is a function of the TLS Bohr radius 
Dt and center separation a^, found using Eq. (fTOl l to be 



-{aT/DTY 



We assume coupling of the TLS to a cold bath in the vac- 
uum state through spontaneous emission, described by a sin- 
gle Lindblad operator L — ^/Va^ ® , where F is the spon- 
taneous emission rate. Measurements of the relaxation time 
of charge qubits in lateral GaAs double dots yielded Ti = 16 
ns,^*" coiTesponding to F w 0.04/ieV~^. 

Transforming Eq. ( |27] | to the interaction picture, with 
p{t) = e'"V(Oe"''"*, we find 



dp 
'dt 



F [2((T_p(T+ — a^a^p — p(T_|-fT_] 



(29) 



with 



,±2 



,±2 



0±2 _ „±2 



± 



+ £L±£2(^r^jQ) + fL_ 



(30) 



In this section we examine the effects of the qubit-TLS cou- 
pling on the performance of a double-dot two-spin qubit. To 
do so we employ a master equation formalism to study the 
dynamics of the coupled qubit-TLS system due to the sponta- 
neous emission of the TLS. Specifically we study the effects 



Here we have defined af — cos flit ± i{uji/fli) sin flit, and 
Ci — {tT/fli) sinftit, where fti = + t^, uji = LOx — a^ 
7 and the upper (lower) sign corresponds to i = 1 (i = 2). 

Inserting a± into Eq. (|29] l we find three separable sets of 
differential equations for the matrix elements of p. The first 
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set is 

^ = -2|a+|4poo + 2ctp22 + 2ci(l - 2c2)Ini(arPo2) 

^ = 2|a+|4poo - 2c\~p22 - 2ci(l - 2cl)lm{a^ po2)m 

^ - - (1 + 2c? - 2c4)/5o2 + 2c2a+2p*2 
+ ia+ci [(3 - 2c?)/5oo + {2cl + l)/522] . 



1 . 
1 . 

1 . 

pPl2 

1 . 

T^P23 



These differential equations are analytically solvable only for 
the case of zero TLS tunneling Ht = 0), where closed ex- 
pressions are also obtained for the matrix elements of the 
original p{t). We include these results in appendix D, as 
they shed light on several features of the general case. For 
nonzero TLS tunneling, equations (l3Tli-(l32li are solved nu- 
merically and the resulting interaction-picture density matrix 
is then transformed back to obtain p{t). In what follows we 
study the effects of a nearby charge fluctuator on a double- 
dot two electron spin qubit by applying equations (|3TTl-([32]i to 
various cases of single qubit rotations. 

B. Single qubit rotations 

In the absence of charge coupling, the time evolution of 
the qubit state is governed by the first term in the Hamilto- 
nian, Eq. ( |28] ). and we obtain a rotation about an axis lying in 
the X — Z plane whose angle with respect to the X axis is 
d = arctan( J/(5/i). When qubit-TLS coupling is introduced 
its deteriorating effects largely divide into the effects of the (3 
and 7 interaction terms. The j3 coupling renormalizes the ex- 
change energy, J ^ J ^ J — 2/3, so that we effectively get a 
rotation about a new axis, whose tilt angle with respect to the 
X axis is I? = arctan[(J — 2/3) /Sh]. 

Figure |5] shows the effects of qubit-TLS coupling on sev- 
eral single qubit rotations. In this figure, as well as for results 
presented in the rest of this paper unless otherwise noted, the 
coupling terms were calculated by averaging over the qubit- 
TLS relative orientations {9,4>,9t,4't) as explained in sec- 
tion lllCI The system's initial state is a singlet for the qubit and 
a localized state for the TLS. Other parameters used are: 
magnetic field inhomogeneity Sh — l/ieV, TLS level splitting 
LOT — 5tT, and TLS Bohr radius and centers half separation 



The equations for {pii, P33, P13} take the same form as those 
above, with ai 02, ci 02- The third set of equations 
reads: 



Dt ~ 5 nm, ay ~ 19.05 nm, respectively. The resulting TLS 
tunnel coupling. It ~ 0.36/ieV, is relatively small thus the 
results can be analyzed in the context of the zero-tunneling 
analytical solution presented in appendix D. 

Figure [5^ shows the singlet probability as a function of 
pulse time for several qubit-TLS distances, R = 50, 100, 200 
nm, at a negative bias point Aa; = 0.382 (below the (1, 1) — 
(0, 2) transition point), where J — 0.15/ieV. In order to re- 
move the effects of the (3 coupling, which are easily corrected, 
we compare the resulting time evolution for each R with a 
rotation about the new tilt-angle The dashed lines corre- 
spond to Ps [t) = 1 — cos^ d SIT? Bt (see appendix D), and 
their deviations from the solid curves correspond to the re- 
maining gate errors due to the combined effects of the 7 cou- 
pling and r. The rotation axis is determined by the (3 cou- 
pling, thus it varies with R. For R — 50, 100, 200 nm we 
find /3 = 0.5, 0.12, 0.03^eV, respectively, and the correspond- 
ing tilt angles are ^ = 48°, 21°, 11°. Figure [Sji shows the 
singlet probabilities vs pulse time for i? = 100 nm at sev- 
eral bias points. Both the exchange and the /3 coupling in- 
crease dramatically as the bias turns positive (see Figure |4| 
thus the corresponding tilt angles for Ax = 0.357, 0.393, 0.4 
are 7? — 6.5°, 52°, 85° respectively, where the latter is close to 
rotation about the Z axis (red-lines and square in figurelS}?). 

Figures [Sj; and |5}l show the pulse error dependence on F, 
corresponding to Figs. and |5}5, respectively. For a self- 
consistent calculation we assume a cubic dependence of F 

with the TLS energy spHtting, F ^ (cj|, + (2iT)^)'^^^, appro- 
priate for free-space spontaneous emission, with the prefactor 
fixed using the data of ref. 56. Thus in figures |5}; and|5}l, F 
variation is accompanied by varying the TLS centers separa- 
tion resulting in a variation in the TLS tunnel-coupling. We 
note that due to the exponential dependence of It on ay large 



I 

= -(2 - c\- c\ - 2ciC2aj"a^)/5oi - ic2 (a^ - 2ciC2a^) po3 + ici{a^ - 2ciC2at)Pi2 + '2c:jclp23 

— mj(c2 + 2ciaj^aJ)poi — {l — c\ + + 2ciC2a^aJ) po3 + 2c?aJ^p^2 + ici{a^ + 2ciC2aJ)p23 

= ia'l{ci + 2c2a^a2)p*Qi + 2c\al'^p*Q^ — [1 + c\ — + 2ciC2a^a2)pi2 + «C2(aJ + 2ciC2ai)p2^ (32) 

= 2a]^^aJ^/5oi + iaj"(2c2aj'aj — ci)po3 — iaj(2ciaj^aj — 02)^^2 ~ (^i + C2 — 2ciC2a]^aJ)p23- 

I 
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FIG. 5: (Color online) (a) Singlet probability vs time for several 
qubit-TLS distances, at bias Ax — 0.382. (b) Singlet probability vs 
time for several bias points at i? = 100 nm. (c) End-of-pulse error 
relative to the /3 corrected rotation about i9-tilted axis vs TLS sponta- 
neous emission rate F, for several R values, at bias Ax = 0.382. (d) 
End-of-pulse error vs V for several bias points at i? = 100 nm. In 
figures (a) and (b) V = 0.04/ieV, the solid lines represent the actual 
time evolution, and dashed lines correspond to pure rotations around 
the i9-tilted axes (see main text). In all the figures pulse duration is 
taken as T = n /5h with magnetic field inhomogeneity 5h = lfj,eW. 



100 nm (7% — 8% error, green circles and dashed lines in 
Figs. |5]3,d), where d is closest to 45° and the F dependence 
is most pronounced. At negative bias both d is close to 0° 
(corresponding to rotation about the X axis) and 7 is very 
small, leading to a small error (blue dotted line in Fig. |5}J). 
At positive bias 1} is fairly close to 90°, but at the same time 
7 is considerably larger, thus the subleading term contributes 
appreciably and the dependence in F is less pronounced (red 
solid line in Fig. |5}l). For yet higher bias, when d — > 90° 
(corresponding to rotation about the Z axis) the smallness of 
the cos^ d in Eq. jP-lj overtakes the increasing value of 7 and 
the errors become extremely small. 

We stress that our results were obtained using a simple 
model for the coupling of the nearby fluctuator with the vac- 
uum, given in terms of amplitude damping. Further investiga- 
tions are required to determine spin qubit dephasing and gate 
errors when other forms of coupling of the fluctuator to the 
charge environment are considered, including specific charge 
noise spectra. 

Notice that single-qubit gates for a two-electron singlet- 
triplet qubit correspond to two-qubit gates for single-electron 
single-spin qubits. Therefore, the results obtained above can 
be directly applied to single-spin qubits as indications of two- 
qubit gate errors. For example, a SWAP gate for single spin 
qubits is done by turning on the exchange splitting J for a 
period of time so that J J/hdt = tt. This coiTesponds to a 
z-rotation for a singlet-triplet qubit when \S + Tq) becomes 
\S — Tq). An error in such a SWAP operation would leave the 
two single-spin qubits with unwanted entanglementi^ 



variations in F amount to a very modest change in ax (8% 
change in ax correspond to three orders-of-magnitude change 
inF). 

Overafl we find that the 7 coupling plays the most impor- 
tant role in determining the pulse errors, while the F coupling 
to the reservoir is less significant. The pronounced dips in the 
pulse eiTor near F ~ 0.1/j,eV, observed in Figs.|5};,d occur due 
to accidental matching between the actual and no-7 rotations 
at the end of the pulse, thus they do not reflect the substantial 
deviations of these rotations that appear throughout the pulse. 
Naively, one would think that for a very small F, the TLS is 
coherent and there is no information loss so that gate errors 
are minimized. While this F dependence is evident in Figures 
|5]: and |5}l, its effect on our results is marginal as compared 
with the gate-eiTors dependence on 7 (notice that the small re- 
duction of ax as F is increased induces a small reduction in 
7 thus the F dependence in these figures is further masked). 
Some insight to the secondary role played by F can be gained 
from the analytical solution for the tx — case given in ap- 
pendix D. 

Inspecting Eq. iP-lj we see that the leading term in the 
expansion of the pulse error in orders of j/B vanishes for 
1) = 0°, 90°, and the subleading term that is still present for 
1) = 0°,is independent of F. Thus the largest pulse eiTors are 
found at Ai = 0.382, i? = 50 nm (14% ^17% error, red 
squares and solid lines in Figs.|5^,c) and at Ax — 0.393, R = 



C. A three-pulse tt rotation about the X axis 

Next we analyze the effects of charge fluctuations on a 
scheme proposed by Hanson and Burkard to produce an effec- 
tive TT-flip about the X axis in the presence of both exchange 
energy and a fixed inhomogeneous magnetic field 6h.^'' This 
setting is desirable since it enables us to perform single-qubit 
operations without relying on a fast control of the interdot tun- 
nel coupling necessary to bring J to zero, an experimentally 
challenging feat. An arbitrary single-qubit rotation can be ob- 
tained with finite exchange energy by applying three succes- 
sive rotations in the X — Z plane. The three-pulse bias cy- 
cle consists of two working points: (i) a negative bias where 
J w Shis small (but need not be zero) and the rotation is about 
a -iJ-tilted axis, (ii) a positive bias in which J ^ dh and the 
rotation is essentially about the Z axis. When Q < -d < tt/4, 
an X rotation at an arbitrary angle ^ can be generated by the 
cycle: 

U,{0^U4x)U^{ip)U4x), (33) 

where Ut){x) and Uz{(p) are rotation matrices about the d- 
tilted and z axes, respectively, and the angles x and (p are 
functions of d and In particular, a 7r-flip about the X axis 
(^ = tt) is obtained when: 

Y = arccosf— tan^ I?); ip — —2 arctan ^ . (34) 
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Without cha rge couplin g, the pulse durations are fixed so 
that x/\/Sh'^ + J|, and = ip/^/SWTTi « (f/Jz, 
where J^, and are the exchange energies at the two bias 
points. Typically most of the duty cycle is spent at the neg- 
ative bias point (i.e., ^ t^). By redesigning the pulse 
sequence to include the renormalized exchange energy , the 
effects of the /3 coupling can be eliminated. The f3 corrected 
pulse sequence is obtained by plugging the new tilt-angle 
■d = arctan[(J^ — 2(3^) /Sh] in Eqs. (l34b . where the d sub- 
script denote the quantity is evaluated in the first (negative) 
bias point. For the parameters we are using, depending on 
the qubit-TLS distance R, the angle-averaged ratio /?/7 can 
be rather small ij3/j = 2.7, 6.6, 16.6 for R 80, 200, 500 
nm, respectively), thus the exchange-renormalization effect 
of 7 is non-negligible. We have found that for short qubit- 
TLS distances (i? < 160 nm, /3/7 < 5), better results are 
obtained by including the 7 coupling in the pulse correction, 
thus this so-called (3' correction is given by using a tilt-angle 
d' = arctan[(J^ - 2(3^ ~ 2j^)/Sh]. 

The above bias cycles were simulated by discretizing bias 
and time steps, so that the actual switching time between the 
working points (which should be nonadiabatic) is taken into 
account. We consider a 7r-flip rotation with the qubit and 
TLS initial states taken, as before, to be a singlet and a lo- 
calized state \(ly[), respectively. Figures |6^-c show the singlet 
probability as a function of cycle time for several qubit-TLS 
distances, R = 500, 200, 80 nm, presenting original (dotted- 
red lines), /3-corrected cycles (dashed-green lines) and /3'- 
corrected cycles (solid-blue lines). The remaining effects of 
the 7 coupling are small (less than 0.1%) for R = 500 nm 
(Fig. |6h), where 7^ = 0.06 neV, but grow rapidly with de- 
creasing qubit-TLS distance. As expected, the two corrections 
deviate only at i? = 80 nm (Fig.|6]:). 

We notice that the three-pulse cycle may fail for one of two 
reasons, depending on the (3 coupling. At large R we have 
\f3z\ ^ Jz 3i the positive bias working point. Then for half of 
the possible system geometries the condition Jz — 2(3z S> 5h 
does not hold and we do not obtain the Z rotation necessary 
to complete the tt flip. For these geometries, however, there 
exists a "sweet spot" where 2(3 — J, thus an effective zero 
exchange can be obtained at this bias, resulting in a single 
pulse X rotation. We shall discuss this case in the next sec- 
tion. At small R (< 30 nm for our chosen parameters with 
5h ~ 1/ieV) a different problem arises due to the large (3 cou- 
pling. Since j3^ ^ the corrected tilt angle d is larger than 
7r/4 for reasonable values of 5h and in fact approaches 7r/2, 
as R decreases. An arbitrary-angle X rotation can only be 
produced when the angle between the axes corresponding to 
the two working points is between 7r/4 and 37r/4, thus, when 
(3 becomes large, the 3-pulse cycle can generate an X rotation 
by a maximum angle of 



arcsin(cot S) f < i9 < ^ 
arcsin(— cot d) 



(35) 



In order to generate a tt flip we need to repeat the three 
pulse cycle, Eq. ( l33b . A^cyc — Floor {^j^j times while re- 
placing the rotation angles given in Eq. ( l34b with Xmax = 
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FIG. 6: (Color online) Singlet probability as a function of time for 
original, /3-corrected, and /^'-corrected pulse cycles, with qubit-TLS 
distance: (a) R = 500 nm, (b) R = 200 nm, and (c) i? = 80 nm. 
(d) TT flip error as a function of qubit-TLS distance. In all figures 
5h^ l^eY, r = 0.04^ieV, Dt = 5 nm, and ar = 19.05 nm. The 
two bias points in the cycle are: Ax^ — 0.357, and Axz — 0.414. 
The kinks evident in these plots are the result of the three pulses 
employed in these gate operations. 



TT. Here, Floor(x) rounds x 



to the nearest integer towards minus infinity. To complete 
the TT-flip, a final cycle is necessary with a rotation angle 
= TT — A^cyc^max, which is generated with^ 



arccos- 



cos t a/ 1 - tan2 d sin^ 4f - sin^ t9 sin^ | 



(/3f = 27r — 2 arctan 



cos^ ^ + cos^ sin^ ^ 



sin xt sin -d 



cos^ ^ 4- cos 2d sin^ ^ 



(36) 



Clearly, as R decreases the f3 coupling becomes larger and 
more cycles are needed to complete the 7r-flip making the 
process inefficient and more susceptible to other decoherence 
mechanisms. Indeed, in this regime, the original premise of 
the bias sequence utilizing two working points with « 6h 
and Jz ^ 6h breaks down and one should consider a dif- 
ferent scheme for performing single qubit rotations. As the 
magnetic field inhomogeneity 6h, used for the qubit rotation 
decreases, this problem will be manifested at larger qubit-TLS 
distances. As far as we are concerned, we push down to the 
small R regime only to demonstrate the increasing effects of 
the qubit-TLS coupling, although it should be noted that our 
results are to be taken with caution when i? < 40 nm due to 
slow convergence of the multipole expansion in this regime. 

Figure |6}l shows the singlet probability at the end of the 
pulse sequence, representing the gate error as a function of 
qubit-TLS distance for original /?— , and corrected cycles. 
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As 7 decreases with increasing R, the difference between the 
two cycles error (green circles and blue triangles) is reduced. 
The non-monotonous behavior of the gate error in the origi- 
nal cycles (red squares) at smaller R appears because in this 
regime the original cycle is so out-of-sync with the actual tilt 
angle -d that its corresponding rotations can accidently bring 
the qubit closer to the triplet state. 

The dependence of the 3-pulse tt flip gate errors on vari- 
ous system parameters is depicted in Figure [T] In order to 
isolate the effects of each parameter on gate errors, we fix 
all other parameters, unlike the calculation presented in fig- 
ure |5] where T,tT,aT were all varied consistently.^'' Figures 
|7^-b show gate error dependence on the size of each of the 
TLS sites (Dt) and TLS half separation (ot), respectively, 
for three qubit-TLS distances (R = 80, 200, 500 nm). While 
D-T is kept fixed in Figure|7j), is scaled with Dt in figure 
|7h, thus the similar behavior of the gate error in both figures 
demonstrates that it is the distance between the TLS centers 
and not their size that impacts the gate error This is consis- 
tent with our identification of the 7 coupling as the source of 
the gate error, since Eqs. dC-St - dC^ show that 7 is propor- 
tional to ttT and does not depend on Dt- At R = 200 nm we 
obtain a gate error of 1% for ot = 175 nm, and a 10% error 
for gt = 964 nm for the /3— corrected cycles. Notice that the 
difference in the remaining error in the two corrected cycles 
gets smaller as R increases. The non-monotonous behavior 
observed in the R — 80 nm case (red lines in Figs.|7h and|7]3) 
are due to higher-order contributions to /? in the multipole ex- 
pansion (PdQ, Pqq) that become significant at shorter ranges, 
and exhibit a complex dependence on Dt (see Eqs. (IC-6b - 
dOT] )). 

Figure |7}; shows the gate error dependence on the TLS 
spontaneous emission rate T. The /3— corrected cycles (solid 
lines) show very little dependence on F, similarly to the X 
rotations studied in the preceding section. Since the cycle in- 
volves d rotations for which F-dependence is present (see Fig- 
ure|5}l and related discussion), we conclude that the F-related 
effects of the first pulse in the cycle are erased by those of 
the third pulse, and the overall gate errors correspond to the 7 
coupling. In contrast, the remaining error in the /3'— corrected 
cycles (dashed lines), which corrects some of the 7-related 
effects, bring out the F dependence. We find the /3' cycles 
error grows linearly with F: Errp/ = aiF for all R values, 
before saturating at a value corresponding to the 7 coupling. 
Saturation is reached at F l/^eV, suggesting that the F dy- 
namics is governed by the TLS tunneling (which is the same 
for all R). For comparison we plot the gate error dependence 
on F for the case of TLS initial state of equal superposition 
+ l^ii)) (dotted lines). This case presents a much 
weaker F dependence since the TLS is set into its ground state, 
limiting the system dynamics. More work is needed to de- 
termine whether these F dependencies are an artifact of our 
simplistic amplitude damping model for the TLS -charge en- 
vironment coupling, or a generic feature characteristic to this 
system. 

Figure |2}1 shows the gate error dependence on the qubit's 
magnetic field inhomogeneity Sh. Below a certain Sh value 
(0.15AieV for i? = 80 nm and O.OT^ieV for R = 500 nm), the 
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FIG. 7: (Color online) tt flip error of /3— corrected (solid lines) and 
/3'— corrected (dashed lines) cycles for three qubit-TLS distances as 
a function of: (a) TLS Bohr radius Dt (ar ~ 4Dt at all points), 
(b) TLS centers half separation ar (Dt = 5 nm for all points), (c) 
TLS spontaneous emission rate F (Dt ~ 5 nm, ar = 20 nm), and 
(d) magnetic field inhomogeneity between the two dots (Dt ~ 5 
nm, ut = 20 nm). The inset in (b) shows TLS tunneling vs. ot. 
Dotted lines in (c) depict the gate error for the /3'— corrected cycle 
with an initial TLS superposition state, (e) Number of cycles needed 
to complete the tt flip as a function of Sh, for /3-corrected cycles, (f) 
Singlet probability vs. time for R — 200 nm, Sh = 0.02/ieV, where 
7 /3— corrected cycles are needed to complete the operation. In plots 
(a)-(c), the magnetic field inhomogeneity is Sh — IfieY and in all 
plots except (c), the TLS spontaneous emission rate is F = O.l^eV. 



f3 corrected tilt angle 1) > tt/4 and more cycles are needed 
to complete the tt flip, as seen in Fig. |7^. In this regime the 
gate error increases for both corrections, obeying a power law 
corresponding to the increasing gate time. Above this thresh- 
old the gate errors of the f3'~ corrected cycles (dashed lines) 
continue to scale with the gate time, whereas those of the 
/3— corrected cycles (solid lines) are largely independent of 
Sh, with a value corresponding to the 7 coupling (and thus to 
R). We find that the scaling of the entire sequence time with 
Sh obeys the power law = aiSh^""^ with fitting parame- 
ters a 1 = 3.23,02 = 1.08 (ai = 3.19, as = 1.06) for i? 80 
nm (R ~ 500 nm), roughly corresponding to a gate time in- 
versely proportional to 5h (slightly larger overhead is needed 
at shorter ranges due to the increased number of cycles). At 
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sufficiently large 5h values error in the Z rotation (2nd step in 
the cycle) is introduced. Since this error grows with 6h its ef- 
fect competes with the 7 related errors, and may thus explain 
the error saturation for the corrected cycles. Apparently, 
the effect of the Z axis error at large 6h is greatly reduced for 
the P'— corrected cycles, where the i9 tilt-angle includes the 
exchange-like contribution from 7, providing a more accurate 
rotation. 

Figure |7f depicts the singlet probability for Sh — 0.02/xeV, 
R = 200 nm, where 7 cycles are needed to complete the tt 
flip, in the /3-corrected cycle (since i5 > 7r/4 in this case, the 
TT flip cannot be completed in the original cycle). We note that 
although the 7 coupling in this case is small (7^ = 1 neV), a 
large gate error of 28.5% is found. This large error is mostly 
due to the small Sh value leading to an overall two-orders-of- 
magnitude increase in the operation time. We stress that for 
short qubit-TLS distances, the large /3 coupling necessitates 
an increasing number of cycles to complete the flip even when 
6h is large. In this case the overall operation times are only 
slightly longer than a single-cycle gate and the resulting large 
gate-errors reflect the large 7 coupling and not the increasing 
number of cycles. 

D. Qubit rotation and dephasing at the J — sweet spot 

The form of the effective qubit exchange energy, J = 
J — 2/3 suggests that certain qubit-TLS geometries can yield 
bias points at which J = 0. At such a bias one can perform 
X rotations in the presence of 6h without resorting to either 
control over qubit tunneling or the three-pulse bias cycles de- 
scribed in the preceding section. In that sense a J = bias 
point is convenient for either performing single qubit rotations 
or doing nothing (idle position). It is thus conceivable that 
in an architecture based on double dot qubit, an additional 
nearby double dot could alleviate the quantum control of sin- 
gle qubit operations. We therefore consider a system consists 
of two double dots, one holding the encoded spin qubit and 
the other assisting in the qubit manipulation, and character- 
ize the effects of the latter on the qubit at the J = working 
point. In general, the f3 coupling can change its sign with R 
but certain geometries yield same sign for /3{R) throughout 
the relevant distance range. For instance, vertically aligned 
double dots (9 = (pr = 0, Or = 90°) yield /3 <0 for any R. 
Here we analyze geometries in which both subsystems lie in 
the same X ~ Y plane. Figure |8] shows four prototypical ex- 
amples, where two configurations (Fig.[8t-b) yield /3{R) > 0, 
and two (Fig.[8};-d) yield (3{R) < for all R. Sweet spots can 
thus be found for the two geometries depicted in Fig. |8^-b. 

To demonstrate quantum control using an ancillary double 
dot, we consider as an example the same axis configuration 
depicted in Fig. [8^, with R = 200 nm. For the second double 
dot (TLS) we take: Dt — 10 nm, ax — 36 nm, = 5 
nm, corresponding to TLS couple-tunneling tr = 0.38/ieV. 
The single-qubit gate is obtained using a single pulse at the 
sweet spot (Aa; — 0.3866 for R = 200 nm). Figure|9^ shows 
the singlet probability vs. time for a 7r-flip rotation with the 
qubit and TLS initial states taken, as before, to be a singlet 




FIG. 8: Several qubit-TLS geometries for 9 — 6t = 90° (qubit and 
TLS lie in the X — Y plane): (a) Same axis (R = xt = £)', (b) 
Parallel axes {R = y,XT ~ x); (c) Perpendicular axes {R = xt = 
y)\ (d) Perpendicular axes (R — x,xt ~ y)- For the geometries 
depicted in figures (a) and (b) the (3 coupling is always positive, while 
for those in figures (c) and (d) it is always negative. 



and respectively. The sweet spot gate operation (solid- 
blue line) presents a much smaller gate error (0.5%) as com- 
pared with the equivalent gate operation obtained using the 
three-pulse /?'— coorrected cycle discussed in the preceding 
section (dashed-green line; gate error 10%). Notice that for 
the rotations shown in this figure the coupling terms are calcu- 
lated with qubit-TLS relative orientations {9, (p, 9t,4't) fixed 
by the specific system geometry and no angular averaging is 
employed. 

Evidently, the J — Q bias point is a convenient idle work- 
ing point. To study qubit dephasing at this point we consider 
Sh ~ for which the qubit computational basis states, S, 
and To are eigenstates of the Hamiltonian, Eq. ( |28l ). The qubit 
initial state is taken as the superposition (IS*) + \To))/\/2, sus- 
ceptible to dephasing (the TLS initial state is still (jy£))- Figure 
|9j' shows the time dependence of the real part of p^j^, the off- 
diagonal element of the qubit's reduced density matrix, corre- 
sponding to the qubit's dephasing (the imaginary part of p^^ 
exhibits a similar time scale). The time scale for dephasing 
is set by the 7 coupling, whereas the role played by the TLS- 
environment coupling T is secondary. For F 7 (dashed- 
green line) and F ^ 7 (dotted-red line), the 7 coupling, which 
entangles the qubit and TLS produces coherent phase oscilla- 
tions. Only when 7 and F are comparable in magnitude damp- 
ing of the qubit phase oscillations occurs (solid-blue line in 
Fig. These results are confirmed by the analytical solu- 
tion of the master equation for the case of no tunneling be- 
tween the TLS centers, given in appendix D. We stress that 
the TLS leads to fluctuations in the energy splitting between 
singlet and triplet states. The randomness in these fluctua- 
tions leading to qubit dephasing comes from the relaxation of 
the TLS to its reservoir rather than from direct interaction be- 
tween the qubit and the reservoir 

It should be emphasized that the behavior described above 
is largely influenced by the choice of parameters, which fixes 
the location of the sweet spot. In particular, since the existence 
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FIG. 9: (Color online) Quantum control using a second double dot in 
same axis configuration (Fig. [8^) with R = 200 nm (a) A vr rotation 
about the X axis with 5h = IfieW, obtained by a single pulse at 
the ,7 = sweet spot, Aa:: — 0.3866 (solid-blue line) and by using 
a three-pulse cycle (dashed-green line). The two bias points in the 
three-pulse cycles are: Aa;^ = 0.357, and Ax^ — 0.414. Here, 
r = 0.04/^eV. (b) Real part of the off-diagonal element of the qubit's 
reduced density matrix vs. time at the sweet spot, for several values 
of r, with Sh = 0. For this figure we used Dt ~ 10 nm, ar ~ 36 
nm corresponding to TLS tunneling rate of tx = 0.38/ieV. 



of a sweet spot requires sufficiently large f3, in order to min- 
imize gate errors and dephasing, one should look for system 
configurations that increase the /3/7 ratio. For the geometry 
considered here we find that gate errors and dephasing exhibit 
a non-monotonous dependence on R, due to the fact that the 
bias point at which J = becomes more negative as R de- 
creases (cl. Fig.|4]l. Since 7 increases for shorter R, but also 
decreases for more negative biases, these two opposing factors 
combine to give a non-monotonous behavior of the gate eiTor. 
While the above discussion has demonstrated the existence of 
sweet spots that may aid in two-spin qubit manipulation, more 
work is needed in order to identify favorable working points, 
so that the proposed architecture becomes viable. 



IV. CONCLUDING REMARKS 

In this paper we studied the effects of a nearby charge fluc- 
tuator (modeled as a TLS) on a qubit encoded by the two- 
electron spin states in a gate-defined double dot. We presented 
a quantitative analysis of the Coulomb coupling between the 
qubit orbital states and the TLS by means of a multipole ex- 
pansion. The resulting interaction terms were divided into the 



/? coupling (/^ X cr^) that renormalizes the qubit exchange 
energy and the 7 coupling (cr^ x f ^) that entangles the qubit 
and the TLS, and can therefore mediate decoherence effects 
to the spin qubit due to the charge environment. We find 
that 7 is generally smaller than (3, and is comparable to the 
exchange energy at i? = 100 nm for the chosen system pa- 
rameters. The ratio p/'j is proportional to both the qubit- 
TLS distance R, and TLS centers separation ax- We note 
that the (3 /j values stated in section IIIICI were obtained by 
averaging over all qubit-TLS orientations, and certain geome- 
tries will yield considerably different values. Both couplings 
depend strongly on R and on the qubit bias. In particular, 
the couplings enhance considerably around the (1, 1) — (0, 2) 
avoided crossing point, where the exchange energy bias de- 
pendence is strongest. Thus, we confirm as expected,-*^ that 
the spin qubit is most susceptible to charge-noise-induced de- 
coherence when operated at positive detuning, at and above 
the anticrossing point. 

As a first application of these results we employed a mas- 
ter equation formalism to study the spin decoherence effects 
due to the charge environment that are mediated by the nearby 
TLS. We derived a set of differential equations for the density 
matrix describing the qubit-TLS system and solved it for var- 
ious scenarios evaluating gate errors during single-qubit op- 
erations, and qubit dephasing times. We find competing de- 
pendence of the gate eiTors on qubit-TLS distance and orien- 
tations, leading to a non-monotonous behavior. For a single- 
pulse rotation, gate errors are found to be largest (5 — 10% for 
r = 0.04/ieV) when qubit rotation is performed about an axis 
midway between the X and Z axes. For positive biases, the 
(3 coupling (at i? < 100 nm) is large enough to produce a Z 
rotation thus, although 7 is large in this regime, gate errors are 
very small. These considerations demonstrate the need for a 
careful analysis of the system at hand in order to optimize its 
performance. 

Analysis of a three-pulse cycle recently proposed to gen- 
erate an X rotation in the presence of finite exchange has 
indicated the implications of a nearby charge fluctuator on 
the feasibility of this scheme. We find that the gate er- 
rors for a TT-flip grow fast with 7, i.e., with decreasing 
qubit-TLS distance. For R = 500, 200, 80 nm the gate er- 
rors are 0.04%, 3.9%, 66% for the /3— corrected cycle, and 
0.02%, 1.6%, 12% for the /?'-corrected cycle, respectively, 
corresponding to 7 = 0.06, 0.88, 13.9 neV at the d working 
point. Moreover for i? < 30 nm, with our parameter choice 
(dh = 1/ieV), the increase in the (3 coupling produces a -d- 
tilt rotation axis that approaches the Z axis, thus more than 
one cycle are needed to complete the tt flip and the scheme 
becomes inefficient. The dependence of cycle errors on other 
system's parameters was also discussed. 

Finally we identified certain qubit-TLS geometries for 
which a convenient working point exists such that the effective 
exchange energy vanishes. We analyze the qubit performance 
at this so-called sweet-spot in terms of gate eiTors and dephas- 
ing {Sh = 0, idle position), focusing on a system geometry in 
which the qubit and the TLS (provided here by an ancillary 
double dot) lie on the same axis. Our results suggest a pos- 
sible qubit design that incorporates a double-dot qubit and an 
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ancillary double-dot that serves to ease qubit manipulation by 
enabling working at the sweet spot. We find that gate errors 
and dephasing are sensitive to the system configuration, and 
in particular to the bias position of the sweet spot, thus a care- 
ful analysis is imperative to optimize qubit performance small 
gate errors. While this design entails an added complexity in 
calibration and initialization of the qubit, it provides an acces- 
sible and convenient working point for an idle position and 
single-qubit rotations. 

In this paper we considered an indirect spin decoherence 
channel where a single TLS's relaxation to a reservoir is 
mediated to the qubit via their Coulomb-assisted entangle- 
ment. Within this simple model, the role played by the TLS- 
reservoir coupling F is secondary to that of 7, and we find 
that pure dephasing appears only when the two are compa- 
rable in magnitude. Our study provides an initial step in the 
quantitative evaluation of the effects of charge environment 
on spin qubits, by microscopically calculating the qubit-TLS 
coupling. To make further progress, one may apply techniques 
developed in studies of superconducting Josephson qubits in 
the context of spin qubits. Notwithstanding the details pertain- 
ing to superconductor devices, fluctuating TLSs weakly cou- 
pled to the qubit were shown to produce both 1// and Ohmic 
noise, inducing qubit relaxation and dephasing4i^Si^ Build- 
ing on these ideas, one can evaluate decoherence effects for 
specific charge noise spectra mediated by the TLS. In addi- 
tion one can consider also direct coupling of the qubit to the 
reservoir and evaluate whether the mediating TLS introduces 
a distinctively new decoherence channel. 

To make contact with current experimental setups of gated 
QDs, one should identify possible candidates for charge fluc- 
tuators, in terms of their physical location and states. Two 
such mechanisms may be electrons jumping between two 
traps or between a localized state and a metallic gate 4^1^ We 
note that in order to produce a random telegraph noise, the 
energy splitting of the TLS should be smaller than kbT so 
that the switching rates for transitions between the two states 
1 — > 2 and 2^1 are comparable. Further measurements of 
the noise spectrum and its temperature dependence, similarly 
to those carried in superconducting devices^^ will be instru- 
mental in pointing at the correct mechanism. 



We find that 
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Here the left (right) pseudospin corresponds to the TLS (qubit) 
state. The concurrence is given by C(V') — \'J2i '^1 \ ^'^'^ Pi'O" 
vides an entanglement measure between the two subsystems. 
Taking the initial states of the two subsystems as a (normal- 
ized) superposition of their basis states: 

= aL\L) + aB.\R) 
l^QiO)) = as\S) + arlT) 

we find the expansion coefficients as: 

ai{t) = asaL(Co-Ci-C2 + C3) + aTai?,(Co + Ci + C2 + C3) 
a2{t) = wsaL(Co-Ci-C2 + C3)-«aTafl(Co + Ci + C2 + C3) 
asit) = -«asaL(Co-Ci + C2-C3)-«aTafl(Co + Ci-C2-C3) 
ai{t) = -a5aL(Co-Ci + C2-C3) + aTafl(Co + Ci-C2-C3) 



and the concurrence takes the simple form: 
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Appendix A: qubit-TLS entanglement due to charge coupling 

Here we evaluate the concurrence of the qubit-TLS system 
state evolving under the interaction Hamiltonian, Eq. ([T]): 



\m) 



-■«.-*|^r(t = 0))® |V^Q(t = 0)> (A-l) 



I leaLfli^asaT (C0C3 - (1(2)! 



SaLttRasarl sin(27i)|. 



(A-3) 



Thus, the qubit and the TLS are entangled only in the presence 
of the 7 coupling and when their individual states are in a 
superposition of their basis states. 



Appendix B: Orbital Hamiltonian for a biased double-dot 

In this appendix we provide details of the calculation of the 
orbital Hamiltonian, Eq. Q within the Hund-MuUiken frame- 
work. This work extends the results of Burkard and Loss^^ to 
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a biased dot configuration. Taking the X-axis (Z-axis) along 
the electric (magnetic) field, the two-electron orbital Hamilto- 
nian is given by 



(B-1) 



where C(ri,r2) = e^/Kjri — is the Coulomb interaction 
between the two electrons, and the single-particle Hamilto- 
nian is 
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The double-dot confinement potential is modeled using a 
quartic potential in the x — y plane and a finite potential in 
the z direction 



yxy(x,y) = 
V.{z) = 
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(B-3) 



where we consider a much stronger confinement in the z di- 
rection, appropriate for typical gate-defines QD structures. 
This enables us to perform separation of variables in the lateral 
and z directions, and to approximate the Coulomb interactions 
using 2-D integrals. The matrix elements of the orbital Hamil- 
tonian are found by adding and subtracting the harmonic po- 
tentials centered at ±a,^^ thus we have 

Horb = /i-a(ri) + /ia(r2) + M^-(ri) + W+ir2) + C 
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Using the orthonormalized single-particle orbitals, ^±a = 
J\f{(j>±a — g(t>^a) we find the single-particle energies and tun- 
nelings in hw^) units: 



6f ±6? 



with 
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Here h is the magnetic compression factor, d, Ax are the dot 
half-separation and electric-field-induced orbital shift in units 
of Bohr radius, respectively, and s is the wavefunction over- 
lap, given in the main text. Notice that the zero-bias results e", 
and coincide with ref. 27- In the basis of the three singlet 
and separated triplet states, the matrix elements of the orbital 
Hamiltonian Eq. Q are found as 
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where the upper (lower) sign in e^[^ corresponds to the sepa- 
rated singlet (triplet) energy. The two-electron Coulomb ma- 
trix elements are 
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(B-8) 



and their closed-form expressions can be found in ref.|2 



Appendix C: qubit-TLS Coulomb coupling terms 



The qubit-TLS Coulomb coupling terms in Eqs. (|24]|-(|26]| 
were calculated to quadrupole-quadruploe order. To 2nd order 
in the qubit orbital overlap, s, they are given by the following 
expressions 
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Appendix D: Analytic solution for the qubit-TLS-reservoir 
master equation for the case of zero TLS tunnehng 

In this appendix we present an analytic solution to the mas- 
ter equation ( |29] l for the case of zero TLS tunneling, <t = 0. 
The a± operators in the interaction picture, Eq. (|30] | reduce in 
this case to 

a± = e±2*("^^") [cos27i(a^ (E) I^) 

± isin27t(cr]^«)crf)] . (D-1) 



The differential equations for the system density matrix as- 
sume a simple form: 
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and are easily solved, with the appropriate initial conditions. 
In order to transform back to p{t), we calculate e^*^*, where 
the system Hamiltonian is given in Eq. ( |28] ) with Bq = 
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0, J), and Bt = (0,0, wt)- The resulting matrix is block diagonal: 
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where we have defined B 



± _ 



-Br 



/3 ± 7, and ri± = 



Bt^. The qubit dynamics can then be examined by 
letting the system evolve under the above Hamiltonian, and 
tracing out the TLS subsystem, thereby obtaining the qubit 
reduced density matrix. 

Eqs. (ID-2I) -( ID-3I ) are used to evaluate the gate error 
of a single-qubit rotation due charge coupling with TLS- 
environment amplitude damping. We perform a rotation about 
the i?-tilted axis, where d ~ arctan( J/(5/i), thus the qubit 
pseudofield vector lies in the X — Z plane. Taking the initial 
qubit state as a singlet, and the TLS initial state as an equal 
superposition of L/R states: 

'/'Q(i-0)= J) ; VT(t = 0) = -^(^}^ (D-4) 

the singlet probability as a function of time is found to be 
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Without the qubit- TLS coupling, the singlet probability at the 
end of a 7r-rotation about the ?9-tilted axis (T^ = ti/2Bq) is 



Inspecting Eq. ( |D-7t we find that the leading term in the gate 
error vanishes for rotations about the X axis {d — 0) or about 
the Z axis {d — 7r/2). Thus the error we expect to observe 
in a TT flip operation will result, predominantly, from the 2nd 
order term, which is independent of F, the latter entering only 
odd order terms in the 7 /B expansion. The F decay contri- 
bution will thus be effective only for rotations about a 6 tilted 
axis. It should be noted that the convergence of Eq. ( ID-7b de- 
pends on the ratio /3/7, thus, as the TLS size increases, 
j3 /7 becomes smaller, and more terms are needed in the 7/-B 
expansion. 



Next we use Eqs. ( ID-2l i-( ID-3b to evaluate qubit dephasing 
at the sweet spot J ~ Q (see section HIlDl i. In order to have 
an idle working point we take 5h = so that B = 0. We 
consider an equal superposition of singlet and triplet (L and 
R) states as our initial qubit (TLS) state. In this case ri± — ^ 
and the equations are trivial with the following solution for 
the off diagonal element of the qubit reduced density matrix: 
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To evaluate the pulse error due to the qubit- TLS coupling we 
design our new 7r-pulse to match the renormalized qubit ex- 
change energy: J J ~ J — 2(3. Thus we have an effective 
pseudospin field B = \{5h, 0, J — 2/3), which redefines the 

tilting angle to i9 = {J — 2/3)/_B, and the corresponding the 
TT -pulse time to — ti/2B. With the effect of the (3 cou- 
pling removed the remaining pulse error is due to 7, and F. 
To second order in 7/i? Eq. (ID-5l l yields 



Ps{T-n) = l- cos^ d + sin 2-d cos 

^2 



cos^ d 



1 



1 - e-^i^/^ 
2 



4 - — 1 sin^ d 
4 



1_ 
B 

(D-7) 



The asymptotic behavior of this equation is: 
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Thus, for our simple amplitude-damping model we find that 
the TLS -environment coupling F has an impact on the qubit 
dephasing only when F « 7. This behavior is demonstrated 
in Figure |9t). 
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